Thermal Spin- Transfer Torques in Magnetoelectronic Devices 
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We predict that the magnetization direction of a ferromagnet can be reversed by the spin-transfer 
torque accompanying spin-polarized thermoelectric heat currents. We illustrate the concept by 
applying a finite-element theory of thermoelectric transport in disordered magnetoelectronic circuits 
and devices to metallic spin valves. When thermalization is not complete, a spin heat accumulation 
vector is found in the normal metal spacer, i.e., a directional imbalance in the temperature of 
majority and minority spins. 

PACS numbers: 72.15.Jf, 75.60. Jk,85.75.-d,75.30.Sg 



Spintronics seeks to exploit the interplay of conduc- 
tion electron charge and spin transport in nanostructures. 
The modulation of an electric current in a spin valve, 
i.e. a normal- metal spacer sandwiched between two fer- 
romagnets, by the relative magnetization directions is the 
essence of the celebrated giant magnetoresistance (GMR) 
effect. Magnetization reversal by the current-induced 
spin-transfer torque in spin valves or tunnel junction 
] has already been applied in memory devices as a low- 
power alternative to Oersted-field magnetization switch- 
mg|. 

Increasing data storage density and access rate is a 
continuing challenge for the magnetic recording industry. 
The relatively high current densities and voltages that are 
required to operate magnetic random access memories 
give rise to heating effects that will complicate modeling 
and deteriorate device stability and lifetime, making it 
difficult to reduce device sizes. Controlled heating can, 
however, also be beneficial: recording by thermally as- 
sisted reversal of magnetization via short laser pulses [3| 
or by Joule heating at highly resistive thin layers Q is a 
possible solution for the next generation of high-density 
non- volatile data storage. 

Johnson and Silsbee [BJ and Wegrowe @ analyzed the 
thermodynamics of transport in collinear ferromagnetic- 
normal metal heterostructures in the diffuse regime. The 
measured magnetothermoelectric power and Peltier ef- 
fect of multilayered nanowires has been described in 
terms of spin-dependent and spin-flip processes in the 
bulk layers by Gravier et al. [7j. A large Peltier cool- 
ing effect in transition metal nanopillars has been mea- 
sured by Fukushima et al. [8]. Nonlinear thermoelec- 
tric transport in non-collinear magnetic tunnel junctions 
has been studied numerically in a tight-binding approx- 
imation 0. Tsyplyatyev et al. invoked thermally ex- 
cited spin currents [1 01 ] to explain thermomagnetic ef- 
fects in metals with embedded ferromagnetic clusters 
llj . In this Letter we report another example of "spin- 



caloritronics" , viz. a strong coupling of thermoelectric 
spin and charge transport with the magnetization dy- 
namics in nanoscale magnetic structures. We establish the 
existence of a thermally induced torque on the magne- 
tization at an interface between a normal metal and a 
ferromagnet. The conditions that should be met in or- 
der to observe this phenomenon experimentally are dis- 
cussed in detail for disordered transition metal-based 
ferromagnet | normal metal | ferromagnet spin valves which 
opens the possibility to switch magnetizations by e.g. 
pulsed laser heating. 

In bulk metallic systems, electron transport is well de- 
scribed by semiclassical diffusion theory [12J. However, 
atomically sharp interfaces should be treated using quan- 
tum mechanical scattering matrices In mesoscopic 
systems such as quantum point contacts scattering 
theory is a well established framework for understanding 
thermoelectric transport (TEl ]. Here we treat magnetic 
nanostructures in terms of electronic distribution func- 
tions in "bulk" layers that are connected with boundary 
conditions at interfaces in terms of microscopic scatter- 
ing matrices us ing an extension of a finite-element (or 



circuit) theory [16|, [17|, [18( to treat charge, spin and en- 
ergy currents on an equal footing. Interface scattering 
is parameterized by a few material-specific conductances 
that are accessible to first-principles calculations. 

We start by partitioning a conducting structure into 
discrete low-resistance nodes connected by resistive ele- 
ments. Ferromagnetic (F) or normal metal (N) nodes are 
characterized by 2 x 2 distribution matrices in spin space 
that can be expanded into a scalar and a vector compo- 
nent / F ( N ) = / C F(N) 1 + & ■ s F ( N )/ s F(N) . The unit vector of 
the spin quantization axis s F is parallel to the magneti- 
zation of the ferromagnet, whereas s N can point in any 
direction. An imbalance between the distribution func- 
tions at two neighboring nodes induces a non-equilibrium 
current. In linear response, the 2x2 spectral current in 
spin space across a ferromagnet-normal metal junction 
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at energy e in the absence of spin-flip and inelastic scat- 
tering is given by Ohm's law 17 1 



«N|F (e) 



Ct/3 



G^(e)u Q / F ( e )-/ N ( e ) 



(1) 



where u^C-D = (1 ± «r-m)/2 are projection matrices in 
which 1 is the 2x2 unit matrix, the unit vector m de- 
notes magnetization direction in the ferromagnet, and <x 
is the vector of Pauli matrices. The conductance ten- 
sor elements read G a " = (e 2 /h) £ nm [<5 TO „ - CJr^ m )*] 
in terms of the energy-dependent reflection coefficients 
r" m (e) for majority and minority spins at the N|F in- 
terface. The total charge-spin and heat matrix currents 

are defined as I = J de ?(e) and eQ — J de(e — fJ-)i(e) = 
eP — fj,I, respectively, where fi is the equilibrium chemi- 
cal potential and I e the energy current. The charge and 
spin electric currents I c and I s are the scalar and vector 
components of the matrix current / = (J c l + a ■ I s )/2. 

Analogously, Q = (Q c l + <x ■ Q s )/2. 

When inelastic scattering in a given node is weak, the 
concept of a local temperature is not applicable and the 
distribution function has to be determined as a function 



of energy [19( , as will be discussed in a future publication. 
Here we assume either that the applied voltage is much 
smaller than the temperature or that there is sufficient 
inelastic scattering so that /r(j.) = fc^fs m ay be pa- 
rameterized by Fermi-Dirac distribution functions with 
spin-dependent chemical potentials MfQ) = //— eVj-Q) and 
temperatures T-m) that are weakly perturbed from their 
values at equilibrium (/i,T). When conductances do not 
vary too rapidly in an energy interval ksT around the 
Fermi level, Sommerfeld's expansion of the distribution 
functions up to order (kuT I ii) 2 may be invoked, where 
kg is Boltzmann's constant [20(. Defining charge and 
spin temperatures T c = (T* + T*)/2 and T s = T T - T l 
(similarly for [i c and [i s , V c and V s ) , we also require T s <C 
2T C in the following. The Sommerfeld expansion leads to 
integrals of the form J de(e — fi) d f s (e) that for d — 0, 1, 2 
read fi s , (ir 2 k 2 B / 3)TT S and (ir 2 k 2 B /3)T 2 /j, s , respectively. 
The same integrals over the function f F — result in 
similar expressions by fj, s — > fi F ~ A*c an d T B — > Tjf — T^ . 
The spin and heat currents through an N|F interface are 
spanned by longitudinal components polarized parallel 
to m (Is = m • I s and Ql = m • Q s ) and transverse 
contributions if = I s — jjm and Qf = Q s — (jlm. 
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The matrix that relates the particle, heat, and spin currents is equivalent to those found in the literature [B|, @] 
when m II s: 
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where G = G^ +G^ is the total conductance, S = —eLoTd e In G\ CF is the thermopower (Mott's law), both at the Fermi 
energy (ep) and L — (7rfcs/e) 2 /3 = 2.45 x 10~ 8 W£!K~ 2 is the Lorenz number. P — (G T — G J -) /G is the polarization 
of the conductance with |P| < 1 and P' = d e [G' — G^) \ ep / d t G\ tF is the polarization of its energy derivative at the 
Fermi energy. In contrast to P, \P'\ is not bounded and P'S can be very large when a van Hove singularity is close to 
the Fermi energy for one spin direction. P$ = (S^ — S±) / (S-\ + Si) = (P' — P) / (1 — P'P) is the spin polarization 
of the thermopower. We focus here on the transverse spin currents: 
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which are parameterized by the (spin-)mixing conduc- 
tance G**- and its energy derivative Gj^ = d e G'^\ 6F . 
We dis reg ard the imaginary part of the mixing conduc- 
tance jl 81 ] and its energy derivative. In analogy with 
the dimensionless mixing conductance rj = 2ReG^/G we 
also introduce a dimensionless "mixing thermopower" as 
if = 2ReGj^/Cc- Both transverse spin currents I, and 
arc absorbed by the ferromagnet and transferred as 
a torque on the magnetization order parameter. 

We extend the methodology used to calculate bare 



interface conductances at the Fermi energy 2l| to ob- 
tain its energy dependence G a ^(e). A finite drift 
is taken into account by replacing G Q ^(e) _1 with 
G^(e)- 1 - (h/2e 2 ) (N N (e)' 1 + S a pN aF (e)- 1 ) [S (H, 
where N aF (e) is the number of propagating modes of 
spin a at energy e in F. To determine the thermopower 
lnG a ^(e) is differentiated numerically. The results for S, 
P' , Ps, T) and r/' are listed for a number of important 
interfaces in Table I. Note that the spin polarization of 
the thermopower in bulk magnets, believed to be domi- 
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|(nV/K 2 ) p'(%) p s {%) n v ' 

Cu/Co(001) -13 72 -8 0.50 0.03 

Cu/Co(001)* -34 89 43 0.49 0.06 




Cu/Co(110) 
Cu/Co(110)* 


-10 
-13 


6 
85 


-66 
45 


0.67 
0.63 


-0.32 
0.07 


Cu/Co(lll) 
Cu/Co(lll)* 


-15 
-15 


56 
77 


-6 
17 


0.53 
0.64 


0.13 
0.13 


Cr/Au(001) 
Cr/Au(001)* 
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Cr/Fe(001) 
Cr/Fe(001)* 


22 
7 


-40 
-190 


48 
-9500 


4.23 
3.25 


-4.27 
-0.48 


Cr/Co(001) 
Cr/Co(001)* 


62 
23 


-111 

-95 


-160 
-92 


3.03 
2.92 


-2.86 
-0.86 



TABLE I: Thermoelectric interface parameters calculated at 
the Fermi energy for a number of almost lattice-matched in- 
terfaces including a drift correction [l3l , [Isl ]. The star * in- 
dicates a dirty interface modeled in a 10x10 lateral supercell 
with two layers of 50%-50% alloy. 



FIG. 1: Schematic of a noncollinear F(mi) \N\F(m.2) spin- 
valve. Spin-dependent resistive elements separate the ferro- 
magnetic reservoirs and the normal metal node. A temper- 
ature bias induces a spin accumulation in the form of heat 
and angular momentum imbalance, the interplay of which is 
governed by inelastic scattering. 



nated by electron- magnon spin-flip scattering 0, [22[ , has 
a different origin. 

The temperature T c , voltage V c , particle spin- 
accumulation sV s and temperature spin-accumulation 
sT s of a given node are governed by Kirchhoff rules. 
Charge and angular momentum conservation implies that 
the sum of all charge and spin currents flowing into a 
given node vanishes. Since thermal transport in metals is 
dominated by the conduction electrons [23j we disregard 
the phonon contribution to the energy currents. Elec- 
trons experience inelastic electron-electron and electron- 
phonon collisions. T s decays by the energy exchange be- 
tween the electrons. The spin accumulation V s is dissi- 
pated to the lattice by spin-flip scattering which can be 
very weak in selected metals and is disregarded here for 
simplicity. We distinguish two different regimes by com- 
paring the dwell time = e 2 T>/ (AttG) , where V is the 
density of states, with te, the energy relaxation time: 
electrons are completely thermalized when te <C tj, 
but effectively non-interacting in the opposite regime 
te Td- The electron dwell time in metallic nanopil- 
lars with a spacer thickness of 10 nm can be estimated 
to be ~ 100 fs. At low temperatures this can be much 
shorter than either electron-electron or electron-phonon 
scattering times 2J| and the spin temperature difference 
or spin heat accumulation becomes an important param- 
eter. For elevated temperatures inelastic scattering is 
more effective and we adopt a complete thermalization 
model. 

We illustrate the theory for symmetric 
Fi(mi)|N|Ffl(m 2 ) spin valves (see Fig. [T]) consist- 
ing of two ferromagnetic reservoirs separated by a 
normal metal node via two resistive contacts with vari- 
able magnetization directions. We calculate the electric 
particle and heat currents and the spin-transfer torques 
for a voltage bias AV = Vr — Vl and temperature 



bias AT = Tr — Tl in the thermalized as well as 
non-interacting regimes. 

In the thermalized regime the spin heat accumulation 
jiN(f) _ q j n ^ e s -j- ea( jy s t a te, lattice and electron tem- 
perature are the same and energy is conserved in the 
node's electronic system. We find for the total electric 
current as a function of the angle 9 between the two mag- 
netizations 



G 



PG tan 2 9/2 
~2~T] + tan 2 9/2 



(AV + SAT)- - ~ ba ^_"t*, n (PAV + P'SAT) . 

(4) 

The angular magnetoresistance for AT = 0, measured by 
Urazhdin et al. [25J, is well described by circuit theory 



26]. When the current bias vanishes, a temperature bias 
AT induces an angular magnetothcrmopower Ay that 
depends on both P and P' 



-AV 



sat) , 



n 



(1 — PP') tan 2 9/2 



(1 - P 2 )tan 2 9/2 



V 



(5) 



The angular dependence of the heat current 

Qc 



-n^(AT/ + 5AT)- |AT 



n 



P'G tan 2 9/2 
~2~T] + tan 2 9/2 



(PAV + P'SAT) 



(0) 



where II = ST is the interface Peltier coefficient, strongly 
violates the Wiedemann- Franz law (n w LqTG). A 
non-negative entropy production rate in the Sommerfeld 
approximation requires |5 maa; | = \[L§ ~ 157 /xV / K [j^ . 

The spin-transfer torque exerted on the magnetizations 
by a temperature difference over the spin valves in the 
thermalized electron regime reads (r = tav + Tat) 



G 



2 r](l + cos9) + (1 - cos 9) 



(PAV + P'SAT). (7) 
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We can understand the similarity of the torque induced 
by the voltage and temperature bias as follows. A tem- 
perature difference over the spin valves initially induces 
different temperatures for the spin species in the normal 
metal node. Since we consider here the strongly interact- 
ing regime, such a temperature difference relaxes quickly 
due to collisions that exchange energy between spin sys- 
tems but conserve the total energy. This is possible only 
by generating a spin current and accumulation that sub- 
sequently induces a torque just as the voltage does. 

The dynamics of the magnetic layers is governed by a 
Landau-Lifshitz-Gilbert equation augmented by the spin- 
transfer torque. We use Slonczewski's estimate for the 
critical current [lj that leads to magnetization reversal 
in metallic ferromagnets, replacing PAV with P'SAT. 
A thermoelectric voltage SAT ~ 100/iV corresponds to 
typical switching current densities of 10 7 A cm~ 2 . Assum- 
ing that a laser pulse provides local heating correspond- 
ing to AT ~ 100K, we require P'S ~ 1//V /K, which is 
not an unrealistic value at room temperature (see Table 
I). When the magnetic layers become thicker, the bulk re- 
sistance and thermopower of the layers dominate. The 
series resistor rule S/G ~ J2 i Si/Gi, where G, and Si ac- 
count for both bulk layers and interfaces in a multilayer 
structure, holds for <C 1. Using the bulk parameters 
by Gravier et al. [28| we estimate that the effective ther- 
mopower can be much higher than a /iV /K, implying a 
strongly increased relative efficiency of thermal magneti- 
zation reversal tat/tav for thicker magnetic layers. The 
conditions for thermal spin-torque switching are presum- 
ably more easily met in spin valves based on magnetic 
semiconductors [29$ . 

In the absence of energy relaxation in the normal node 
, a vector spin heat accumulation sT s N develops. When 
S 2 <C Lq, the "non- interacting" thermal spin-transfer 
torque t at reduces to the simple expression: 



t at tat 



TAT 



AV=Q 



(t]'-v)P tan 2 0/2 
rjP' t? + tan 2 0/2' 



_ s 

AT 



r)P' 



'AT 



TAT 



(8) 



(9) 



av=o (V ~ V) sin 6/2 t at 
In conclusion, we presented a circuit theory of thermo- 
electric transport in non-collinear spin-valves. In thinly 
layered structures, transport properties are governed by 
interface conductances and their energy derivatives that 
have been computed from first principles. We predict 
a spin-transfer torque associated with purely thermal 
currents that can be large enough to reverse magneti- 
zations. The concepts of spin heat accumulation and 
spin-mixing thermopower have been introduced to de- 
scribe the thermoelectric transport in different energy 
relaxation regimes. We expect that a temperature gradi- 
ent can excite magnetization dynamics in magnetic tun- 
nel junctions and domain walls in ferromagnetic wires as 
well. 
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